Giant Magnetoresistance in a Chemical Vapor Deposition Graphene Constriction

Magnetic field-driven insulating states in graphene are associated with samples of very high quality. Here, this state is shown to exist in monolayer graphene grown by chemical vapor deposition (CVD) and wet transferred on Al2O3 without encapsulation with hexagonal boron nitride (h-BN) or other specialized fabrication techniques associated with superior devices. Two-terminal measurements are performed at low temperature using a GaAs-based multiplexer. During high-throughput testing, insulating properties are found in a 10 μm long graphene device which is 10 μm wide at one contact with an ≈440 nm wide constriction at the other. The low magnetic field mobility is ≈6000 cm2 V–1 s–1. An energy gap induced by the magnetic field opens at charge neutrality, leading to diverging resistance and current switching on the order of 104 with DC bias voltage at an approximate electric field strength of ≈0.04 V μm–1 at high magnetic field. DC source–drain bias measurements show behavior associated with tunneling through a potential barrier and a transition between direct tunneling at low bias to Fowler-Nordheim tunneling at high bias from which the tunneling region is estimated to be on the order of ≈100 nm. Transport becomes activated with temperature from which the gap size is estimated to be 2.4 to 2.8 meV at B = 10 T. Results suggest that a local electronically high quality region exists within the constriction, which dominates transport at high B, causing the device to become insulating and act as a tunnel junction. The use of wet transfer fabrication techniques of CVD material without encapsulation with h-BN and the combination with multiplexing illustrates the convenience of these scalable and reasonably simple methods to find high quality devices for fundamental physics research and with functional properties.

* sı Supporting Information ABSTRACT: Magnetic field-driven insulating states in graphene are associated with samples of very high quality. Here, this state is shown to exist in monolayer graphene grown by chemical vapor deposition (CVD) and wet transferred on Al 2 O 3 without encapsulation with hexagonal boron nitride (h-BN) or other specialized fabrication techniques associated with superior devices. Two-terminal measurements are performed at low temperature using a GaAs-based multiplexer. During high-throughput testing, insulating properties are found in a 10 μm long graphene device which is 10 μm wide at one contact with an ≈440 nm wide constriction at the other. The low magnetic field mobility is ≈6000 cm 2 V −1 s −1 . An energy gap induced by the magnetic field opens at charge neutrality, leading to diverging resistance and current switching on the order of 10 4 with DC bias voltage at an approximate electric field strength of ≈0.04 V μm −1 at high magnetic field. DC source−drain bias measurements show behavior associated with tunneling through a potential barrier and a transition between direct tunneling at low bias to Fowler-Nordheim tunneling at high bias from which the tunneling region is estimated to be on the order of ≈100 nm. Transport becomes activated with temperature from which the gap size is estimated to be 2.4 to 2.8 meV at B = 10 T. Results suggest that a local electronically high quality region exists within the constriction, which dominates transport at high B, causing the device to become insulating and act as a tunnel junction. The use of wet transfer fabrication techniques of CVD material without encapsulation with h-BN and the combination with multiplexing illustrates the convenience of these scalable and reasonably simple methods to find high quality devices for fundamental physics research and with functional properties. KEYWORDS: multiplexed device arrays, graphene, CVD, giant magnetoresistance, magnetotransport, Fowler-Nordheim tunneling, direct tunneling I nsulating phases in graphene at high magnetic fields are typically observed only in exceptionally high quality devices, characterized by high carrier mobility and low residual carrier density. Demonstrations have been almost exclusively limited to exfoliated graphene flakes 1−6 and, recently, ultrahigh mobility graphene 7,8 grown by chemical vapor deposition (CVD) and encapsulated with hexagonal boron nitride (h-BN). 9 These devices are fabricated from monolayer material using sophisticated methods that limit scalability such as the use of exfoliated material, suspension of the graphene, 10,11 fabrication on h-BN substrates, 12,13 or encapsulation with h-BN. 9,14,15 At a high perpendicular magnetic field, the 2D density of state becomes quantized in Landau levels with energy = ℏ | | E Nv e BN sgn( ) 2 LL F , where N is the Landau level index, v F is the Fermi velocity, ℏ is the reduced Planck constant, and e is the electron charge, and an insulating phase can exist at charge neutrality 3,14 characterized by a resistance that diverges with B as a magnetic fielddependent energy gap opens in the N = 0 Landau level. Here, we find such an insulating state can exist in conventionally wet transferred monolayer CVD graphene devices fabricated on amorphous oxides.
The device is contained within an array of CVD graphene devices nominally patterned as 10 μm-wide squares measured using a multiplexer 16 at cryogenic temperatures. Multiplexing is a means for high-throughput testing and characterization of many devices, 17−19 and the large number tested can lead to the discovery of individual devices with exceptional qualities that can be opportunistically studied. Here, a narrow constriction of ≈440 nm is created at one end of the graphene square, apparently by graphene tearing during fabrication. The resistance at charge neutrality diverges with B showing insulating behavior and an energy gap that increases with B. The conductivity appears to be reasonably temperatureindependent below ≈2 K, and an increase in T shows activated behavior from which the energy gap is estimated to be ∼2.4−2.8 meV at B = 10 T. This behavior suggests the existence of a region of high electronic quality within the narrow channel, which dominates transport at high B. The analysis of charge transport mechanisms using DC source− drain bias measurements and quantum tunneling models 20−24 at low T shows a transition from direct tunneling to Fowler-Nordheim tunneling with increasing DC bias from which the length of the tunneling region is estimated to be on the order of ≈100 nm. The device shows the current switching of order ∼10 4 with DC source−drain bias for an approximate electric field strength of ∼0.04 V μm −1 .
Here, although the narrow constriction is unintentionally created, it is straightforward to create the same geometry deliberately. A second device is fabricated with a similar contact configuration and a 600 nm-wide constriction defined next to one contact by electron−beam lithography. Similar behavior is observed in terms of an insulating state although being impacted by lithographic patterning and doping, potentially indicating the quality of the constriction edge as an important consideration. We speculate that less disorder is present in the unintentional constriction due to its formation as a result of tearing, where edges are never exposed to plasma etching, which can cause greater edge disorder and contamination compared to cleaner alternatives 25,26 thus motivating further studies toward developing clean processes for constriction edges such as h-BN etch masks. 27,28 The observation of the insulating phase at charge neutrality in wet transferred CVD graphene on amorphous oxide suggests the intriguing potential of developing systems to investigate exotic quantum phenomena using conventional and scalable fabrication techniques. A primary application of multiplexing is the statistical analysis of large data sets and identifying trends in data ensembles. 16,29,30 Statistical analysis identifies general behavior, allowing outliers to be selected and studied in detail. Our present work highlights the capability of our multiplexing system for an in depth study of the physics of such channels using numerous techniques including magnetic field, source− drain bias, and temperature from which one can obtain a detailed understanding of the fundamental phenomena. The identification of optimal devices has been a bottleneck in studies of fundamental physics, since multiple devices usually must be measured before a suitable one is found. The multiplexer expedites this process in a much more efficient manner. Our study also highlights the ability of the multiplexer to drive device development, since the inherent inhomogeneity in nanostructure devices leads to these outliers that can have desirable properties. The study of the outliers reveals a particular set of device conditions creating the observed behavior, and the identification of differences from other devices in the array provides information needed to reproduce the behavior deterministically.

RESULTS AND DISCUSSION
The multiplexer is fabricated on a GaAs/AlGaAs heterostructure containing a two-dimensional electron gas (2DEG) with a conducting path from one input to 16 outputs. Individual graphene devices are connected to each output and selected using addressing gates. Multiplexer operation, fabrication, and layout are described in ref 16. The device measured here is different from that in ref 16 where the emphasis was on the versatility of the multiplexer for highthroughput testing and compatibility with diverse nanoelectronic devices via integration of arbitrary nanomaterials including 2D materials (both CVD-grown films and exfoliated flakes) and semiconductor nanowires. 16 The multiplexing platform can also be adapted for devices fabricated within a 2DEG and controlled by surface gates. 19,29−33 Here, the array of graphene devices is patterned in a monolayer CVD graphene film after wet transfer 34 to the multiplexer host chip, as described in the Methods. Scanning electron microscopy images of the measured device are shown in Figure 1a,b. The image shows a narrow constriction next to a contact, and the region around the constriction appears to have reasonably well-defined edges and angles, suggesting the graphene constriction was formed by tearing along clear crystallographic directions. Tearing has been observed to give armchair or zigzag termination, 25 which could indicate the constriction edges are atomically clean and feature minimal dangling bonds and scattering sites. 26 The constriction is ∼440 nm at the narrowest point. Figure 1c shows a schematic cross section through the device (side profile). Source−drain contacts are created prior to graphene transfer, and 95 nm thick Al 2 O 3 is deposited by atomic layer deposition above the back gate. The two-terminal differential conductance is measured in a 3 He cryostat.
Transport Characterization. Figure 1d shows resistance (R) as a function of back gate voltage (V G ) at B = 0 T and T = 0.29 K. The black curve is a fit to the data using 35,36 i k j j j j j j j y where R P is the parasitic resistance, L/W is the length to width ratio, μ is the device mobility, n 0 is the residual carrier density, and n is the back gate-dependent carrier density. The density is given by n = C G (V G − V 0 )/e, where C G is the gate capacitance per unit area C G = ϵ 0 ϵ/t ox , V 0 is the charge neutrality-point voltage, ϵ 0 is the free space permittivity, ϵ is the oxide permittivity, and t ox is the oxide thickness. The mobility and parasitic resistance are estimated separately for electrons and holes, i.e., for V < V 0 (V > V 0 ), μ = μ hole , and R P = R P,hole (μ = μ electron and R P = R P,electron ). For this equation, V 0 is obtained by initially fitting data with density-independent R P and μ, giving 0.52 V. Fitted parameters are R P,hole = 17.4 kΩ, R P,electron = 19.0 kΩ, n 0 = 1.3 × 10 11 cm −2 , μ hole = 5350, and μ electron = 6360 cm 2 V −1 s −1 for L/W ≈ 2.3. L/ W is estimated by dividing the conducting area into squares and analyzing the resulting resistor network to give an equivalent length-to-width ratio. Figure 1e shows the Landau level splitting with the magnetic field 37−39 with dR/dV G plotted as a function of B and V G . Local resistivity maxima corresponding to where ±1 Landau levels are half filled are indicated by dotted white lines, and regions between correspond to quantum Hall plateaus. The gate capacitance is estimated from the gate dependence of quantum Hall states ν = ±2, where ν is the filling factor, using dn/dB = νe/h and dn/dV G = C G /e, giving C G ≈ 60 nF cm −2 and ϵ = 6.4. Landau level splitting also provides an alternative estimate of μ using the approximation of μ ≈ 1/B at B where Landau levels can be resolved. 40 This gives a possible mobility up to μ ≈ 8000 cm 2 V −1 s −1 , as discussed in Figure S3a, and mean free path l = v F τ ≈ 17 nm using scattering time τ ≈ 1/ω c , Magnetoresistance Measurements. Figure 2a shows R as a function of V G at B = 7, 8, 9, and 10 T and T = 0.29 K. The resistance at charge neutrality increases by almost 3 MΩ at B = 10 T. Figure 2b shows R at V G = 0.31 V as a function of T, which diverges with B for the lowest T. The estimated Bdependent resistance of the GaAs multiplexer is subtracted from the data in Figure 2 as described in Figure S1. Data in Figure 2b,c are measured at V G = 0.31 V. Data for Figure 3 are obtained at V G = 0.32 V after reoptimization. Both are taken to represent the high B insulating state at ν = 0 when discussing each figure. Fitting transfer characteristics at B = 0 T using eq 1 estimates a charge neutrality point voltage of 0.52 V. Differences may arise since the low B value is given by fitting methods over a broad range of V G , whereas values at B = 10 T are experimentally determined to be as close as possible to the maximum R. Additionally, differences may conceivably arise since the low B characteristics likely represent an average of the entire graphene area, whereas at high B a small region within the constriction appears to dominate transport characteristics.
The conductance G at V G = 0.31 V is plotted as a function of T on a lin-log scale in Figure 2c for B = 7, 8, 9, and 10 T. This behavior becomes visible for B ∼ 7 T and above. Similar to that in ref 1, G appears to become relatively T independent at the lowest temperatures measured. Solid lines show fits to the data for activated transport using a Fermi−Dirac distribution 15,41 where E a = 0.7, 1.0, 1.6, and 2.4 meV is the magnetic fielddependent energy gap at B = 7, 8, 9, and 10 T, respectively, k B is the Boltzmann constant, γ is the zero T conductance, and β is a pre-exponential factor. Fits at each B are performed independently with E a , γ, and β as fitting parameters. Data are first fit using the Arrhenius equation, which gives energy gaps comparable to the thermal energy k B T at higher T values for which the measurement is performed (up to T = 25 K); therefore, the Fermi−Dirac function is used instead. Fitting using the Arrhenius equation is shown in Figure S2b.
The energy gap is shown as a function of B in Figure 2d. If a linear dependence E a = ΔE Z − Γ is assumed, 5,15 where ΔE Z is Zeeman energy gμ B B, μ B is the Bohr magneton, and Γ is the disorder broadening of Landau levels, parameters Γ = 39.2 K and effective Landég-factor g ≈ 10 are enhanced above the free-electron g-factor g = 2, as seen previously for exfoliated graphene, 2,13 and attributed to exchange interactions. 15, 42 The quantum lifetime 5,43 estimated from Landau level broadening The fitting parameters using eq 2, i.e., energy gap E a , fitted conductance at zero temperature γ, and pre-exponential factor β, are shown in panels (d), (e), and (f), respectively.
is reasonably comparable to scattering lifetimes estimated from mobility using the Drude formalism τ = (hμ/2ev F )(n/π) 1/2 ≈ 105 fs at n = 2 × 10 12 cm −2 . However, the range of data in Figure 2d is too small to give a true indication of the relationship between E a and B. An approximately linear relationship appears in some devices described in the literature 2,13 but not necessarily in others. 10 Reference 14 shows a linear relationship for a device engineered to minimize the Coulomb interaction and a square root dependence consistent with a Coulomb energy relationship in a separate device without this reduction. The estimated gap sizes here are much smaller than the Coulomb energy . Previous studies show either E a to be small 2 compared to E C or that they are similar values. 13 Fitting parameters γ and β are shown in Figure 2e,f, respectively. The γ value describes G in the low T limit and nears zero as B increases. In a scenario of transport dominated by a region of gapped graphene at high B, this may represent hopping transport that diminishes as the energy gap E a increases and the size of the wave function reduces with B, leading to less hopping between sites. The pre-exponential factor β is reasonably similar for B = 8, 9, and 10 T for fits performed independently at each field and is reduced at lower B as the gap becomes small. The β term may be related to the geometry of the gapped region, becoming reasonably constant above a crossover field where the gap opens. The insulating behavior does not occur in other devices in the multiplexed array that do not contain a narrow constriction. The previous observation of insulating states in graphene associated with high electronic quality devices supports a scenario where a high quality region exists within the constriction and dominates the transport characteristics at high B.
DC Source−Drain Bias Measurements. The insulating state is further investigated by applying a DC source−drain bias across the graphene device (V DS ). Figure 3a shows R as a function of V DS at select V G around charge neutrality for T = 0.29 K and B = 10 T. Data representing the high B insulating state at ν = 0 are obtained at V G = 0.32 V; data at other V G are labeled according to Δn = C G (V G − 0.32)/e. Supporting Information describes how V DS across the graphene is extracted from the total source−drain bias applied to the circuit (V B ), accounting for series resistance. The temperature dependence at ν = 0 is shown in Figure 3b. Energy scales associated with DC bias are reasonably similar to E a from Fermi−Dirac fitting, i.e., in Figure 3b, the full width at halfmaximum (fwhm) of ≈3 mV at T = 0.29 and 2 K, accounting for a vertical offset from series resistance (R S ) of ≈219 kΩ. Data in Figure 3a,b are plotted before subtracting R S . There is a small effective gating of the graphene at high V DS due to the source−drain bias, discussed in Figure S4; however, the effect is negligible. Figure 3c shows G as a function of T at select V DS after subtracting R S . For clarity, data are shown for V DS ≤ 0, and data for V DS > 0 are similar. The dashed line shows a fit to V DS = 0 mV data using the Fermi−Dirac distribution. Parameters are as follows: E a = 2.8 meV, α = 0.02 e 2 /h, and β = 1.4 e 2 /h. Differences from the parameters estimated for B = 10 T data in Figure 2c may arise from reoptimizing V G to find the maximum R. Data for finite V DS appear to fall along the V DS = 0 mV fit line at high T and diverge as T is reduced, becoming reasonably T independent for T < 2 K. The gap ΔT between T = 2 K data points for finite V DS and the V DS = 0 mV fit line is shown in the inset converted to energy ΔE = k B ΔT. Data at T = 2 K are chosen since G becomes reasonably T independent below this T. This may represent an equivalent thermal energy or an equivalent reduction in barrier height that would be required at V DS = 0 mV to achieve the same effect. Dashed lines show linear fits with the same gradient (absolute value) for V DS > 0 and V DS < 0 mV, which cross the horizontal axis at | V DS |≈ 0.6 mV, similar to the AC excitation voltage V AC = 400 μV (rms). The relationship appears linear up to bias voltages similar to the activation energy. We do not observe a power law relationship between ΔT and applied power at each V DS , which could be expected if ΔT represented a rise in electron temperature. 44 This suggests the main impact of V DS within this regime is to change the bias across the barrier. However, this is better investigated by analyzing the current−bias voltage characteristics.
I−V Characteristics. Figure 4a shows the DC current I DS as a function of V DS for select V G at T = 0.29 K. The current is obtained by integrating the measured differential conductance G = dI DS /dV B over the total DC bias applied to the circuit V B . The I DS remains near zero over the widest range of |V DS | ≲ 3 mV at ν = 0, an energy scale similar to E a . Figure 4b shows the temperature dependence for the ν = 0 insulator. Above T = 2 K, data become much more linear, as transport becomes thermally activated. The I DS /V DS characteristics are shown on a log-lin scale in Figure 4c,d, for different n at T = 0.29 K and different T at ν = 0, respectively. Positive V DS data are plotted  Figure S5. The current changes by nearly 4 orders of magnitude for ΔV DS from 0 to ≈4 mV at T = 0.29 K and ν = 0. If the barrier length is on the order of ≈100 nm (see later analysis); this corresponds to an electric field strength of ∼0.04 V μm −1 . The narrow constriction connects to the Ti/Au contact on one side and the wider graphene area on the other. The similarity of data for positive and negative V DS , Figures 4  and S5, respectively, suggests that any asymmetry caused by different types of contacts to a tunneling region within the constriction is small. Data at T = 0.29 and 2 K cross at high V DS , likely caused by uncertainties in the series resistance correction. Therefore, only trends at separate T will be considered. In this measurement, R S depends on T and total bias V B applied to the circuit. The V B dependence of R S arises from a gating effect where the bias across the multiplexer changes the Fermi energy in the GaAs 2DEG relative to the addressing gates. The resulting change in resistance manifests as an approximately linear background with DC bias voltage, dependent on T, and is corrected for as described in the Supporting Information. The accumulation of estimates makes comparing data between different T's less reliable when V DS is the highest and the graphene resistance has become small. It is straightforward to remove the R S dependence on DC bias with changes to the multiplexer operation during the experiment by maintaining a fixed potential difference between addressing gates and the 2DEG. 45 This will allow easier characterization of I−V curves in the future.
Tunneling mechanisms are investigated by analysis of the current−bias voltage characteristics using eqs 3 and 4. 21−24 Equation 3 describes direct tunneling through a trapezoidal barrier at low bias where ϕ is the barrier height, m is the effective mass, and d is the length of the tunnel barrier. Plots of ln(I/V 2 ) as a function of ln(1/V) show a linear relationship in this regime. At high bias (V > ϕ), the barrier becomes triangular and transport can be described using Other studies have also observed eV t to be smaller than the estimated barrier height, 21 attributed in part to the approximation of the barrier being simple trapezoidal/ triangular shapes instead of being fully treated.
The T dependence of the ν = 0 state is shown in Figure 4f,h. Results at T = 0.29 and 2 K are consistent with a transition between direct tunneling and FN tunneling. No inflection  Figure 3c. The barrier length is consistent with the estimated Fermi wavelength λ π = = n 4 / 97 nm F at n = n 0 . The barrier length of ≈100 nm is only given as an approximate estimate, since within the fitting range, the data are not yet truly linear, and it may be that V DS is not yet sufficiently high. Further deviations in linearity can occur due to T and image effects 20,21,47 that are not included in eq 4. The equation also assumes a simple triangular barrier shape and is an approximation in the high bias and zero T limit. These factors can contribute to uncertainties in the estimate of d. Additionally, uncertainties in R S at high V DS may also cause deviations from linearity. The transition between Fowler-Nordheim and direct tunneling has not previously been reported in this insulating regime, and to our knowledge, detailed I−V measurements and analysis including temperature dependence have not hitherto been performed. Where I−V sweeps were carried out at a single T, the nonlinear regime indicated an energy scale an order of magnitude larger than the activation energy, suggesting charge transport may occur over several connected insulating areas. 10 In contrast, our results indicate energy scales similar to E a , perhaps suggesting an individual insulating region.
An aspect that distinguishes this work is the use of bottom contacts. These have potential to induce strain in graphene next to the contact edge, creating the possibility for the physical origin of these effects to lie in the combination of narrow constriction and strain, where strain from metal contacts can affect device properties. 48 The high quality region may also indicate cleaner interfaces and low defect density within the narrow constriction, and strict control of the interface quality may be necessary for high yields. We do not observe symmetry breaking and ν = ± 1 states as have been reported in samples with much higher mobilities 5,49 or measured at much higher magnetic fields. 5 Similar to refs 2, 3, and 50, this may be attributed to comparatively higher disorder in our devices or insufficiently high B. For comparison, devices hosting the high B insulating state with very similar mobilities to ours, refs 2 and 3, did not detect ν = ± 1 states for B up to 30 and 18 T, respectively. It is also likely that the amount of residual carrier density or disorder in our device precludes the existence of low B insulating transitions that appear within the global phase diagram for 2D systems in a magnetic field, 51 where increasing density can reduce the critical field of the transition to B = 0 T. 52 Instead, our results appear to be consistent with the phase diagram in ref 3.
The insulating state is often observed as B is increased, 1−4,10,11,53 despite the fact that the ground state at ν = 0 is predicted to be a quantum Hall ferromagnet with counter-propagating edge modes, where each edge section contributes R = h/e 2 and total R is given by the combination of these edge sections. 13,14,54−56 It has recently been shown that quantum Hall ferromagnetism can exist in devices engineered to screen the Coulomb interaction, whereas the insulating ν = 0 state is present in devices where this is not the case. 14 This was attributed to the indirect modification of lattice−scale electron−electron interactions, which determine the ground state and can lead to either ferromagnetic or insulating phases. 57 The observation of the insulating state, previously requiring the electronic quality given by exfoliation and encapsulation methods, in unencapsulated and wet transferred CVD graphene may be pertinent for the development of graphenebased devices to investigate physics phenomena due to the scalability and relative simplicity of these methods. To address this potential, the role of constriction width and possible strain at the contact edges should be investigated, since here, the behavior appears in a device with a narrow constriction and back contacts. Although the constriction is created unintentionally in the initial device, it is straightforward to determinately recreate with standard lithographic techniques. We deliberately create a similar geometry using electron−beam lithography and oxygen reactive ion etching to define a 10 μmwide graphene channel with a 600 nm wide constriction next to one contact electrode in monolayer CVD graphene transferred on a doped Si wafer. The device layout is shown in Figure 5a−c, and fabrication and measurement details are given in the Supporting Information. An ∼50 nm Al 2 O 3 encapsulating layer is added by atomic layer deposition (ALD) after patterning, 58 and all data presented includes contact resistance. Figure 5d shows R as a function of V G at B = 0, 3, 6, 9, and 12 T and T = 1.5 K. Two clear resistance peaks are visible at B = 0 T (bold trace), labeled a and b. The colormap in Figure 5e reveals their evolution into two sets of Landau levels with charge neutrality points V CNP,a = ∼42 V and V CNP,b = ∼18 V for a and b, respectively, marked by the horizontal dashed lines and arrows. When one considers the graphene channel as two resistive components in series with distinct aspect ratios representing the constriction and main channel, there are indications that the behavior of peak a (peak b) is consistent with arising from the constriction (wide region). Specifically, the higher R of peak a at B = 0 T is consistent with a larger L/ W ratio, which occurs in constriction, as supported by fitting the peaks separately using eq 1, shown in the Supporting Information. Additionally, nanopatterning graphene (the use of PMMA and planar etching) has been shown to lead to pdoping, which may explain the positive shift of V CNP for peak a due to the enhanced impact of these doped etched edges on transport through the constriction relative to the wide section. Dual peak behavior is not observed in the unintentional constriction, possibly due to the tearing process resulting in cleaner constriction edges and reduced contamination such that there is no difference in V CNP for graphene in the constriction and elsewhere in the device. We speculate that removing the larger graphene area, i.e., fabricating only narrow graphene channels between micron-spaced source−drain electrodes, may remove the secondary Dirac peak. Figure 5g shows R at ν = 0 for peaks a and b as a function of B at T = 1.5 K. The resistance of peak a rises with B, whereas b plateaus at ∼18.4 kΩ. The T dependence of peaks a and b at B = 12 T is shown in Figure 5h. Peak a shows insulating behavior, whereas peak b is relatively unchanged in comparison. The behavior of peak a is similar to the device with the unintentional constriction in terms of increasing R with B and insulating T dependence. The magnitude of the change in resistance for the insulating behavior is reduced compared to the unintentional constriction, which can be attributed to a combination of the higher base temperature (1.5 K vs 300 mK) and increased disorder and doping in the etched edges. 27,28 It is expected that devices with higher disorder will show a later onset of the insulating state such that R is lower at the same B and T. 1,11,53 For comparison, Figure  2b shows R = ∼860 kΩ at B = 10 T and T = 2 K. We speculate that less disorder is present in the unintentional constriction since it forms as a result of tearing, 25,26 such that the edges are never exposed to plasma etching. However, the trends of R with T and B in both devices are the same, and fitting data from Figure 5h using eq 2 estimates an energy gap E a of ≈0.9 meV at B = 12 T, compared to E a of ∼2.4−2.8 meV for the unintentional constriction.
We note that the qualitative difference in the B and T dependence of peaks a and b suggests behavior that is inherent to each charge neutrality point, and the effect of overlapping Landau levels is small. If this were not the case, both peaks might be expected to show the same behavior, since peak b overlaps a and vice versa. The strong T dependence at peak a in Figure 5f compared to the little discernible change around peak b and for V G < 18 V (i.e., in the region of Landau levels from peak b) also suggests the T dependence around V G = ∼42 V (V CNP,a ) can be attributed to the ν = 0 state of peak a rather than overlapping Landau levels.
These data provide evidence that the insulating state at ν = 0 can be reproduced in deliberately fabricated CVD graphene devices and that the unintentional constriction may be benefiting from higher electronic quality due to the presumably torn edges. This suggests the importance of lithographic patterning processes that minimize edge disorder and contamination, such as using h-BN masks, 27,28 when creating deterministically defined constrictions, and may imply the possibility for systems to investigate exotic physics using relatively simple and scalable fabrication methods with the implementation of such techniques. Our results highlight the use of multiplexing in driving device development, since by studying the device conditions leading to the initial observation of the insulating state, we are able to reproduce it deterministically and identify key factors for its existence. Additionally, the multiplexing platform combined with scalable fabrication such as wet transfer techniques of CVD material has demonstrated a convenient means by which devices can be identified for fundamental physics research or with functional properties.
Future expansion of the multiplexing tool to achieve fourprobe capabilities and incorporate h-BN substrates or h-BN encapsulation presents intriguing possibilities to explore other diverse phenomena, including those in graphene devices on patterned h-BN. 59 The N = 0 Landau level in hybrid graphene systems has also shown intriguing behavior and could be investigated. 60,61 Expansion to other nanomaterials is straightforward since nanowires, exfoliated 2D materials, and CVDgrown 2D materials can be transferred to the multiplexer using compatible fabrication methods. 16 The onset of the high resistance state studied here in graphene has been shown to occur at lower magnetic fields in devices with higher mobility and lower disorder. 1,11,53 The reduction of the magnetic field required may bring CVD-grown graphene devices into the realm of giant magnetoresistance devices and magnetic switches 62 and motivate the development of graphene-based magnetic field detection or magnetic texture mapping sensors with cryogenic readout.

CONCLUSIONS
These results show the existence of a magnetic field-induced insulating state in monolayer CVD graphene with resistance that diverges with B at charge neutrality. This state is usually associated with very high quality devices made using exfoliated graphene flakes and suspended, encapsulated, or fabricated on extremely smooth substrates such as hexagonal boron nitride. The device studied here uses CVD-grown graphene that is unencapsulated and fabricated on an Al 2 O 3 substrate created by atomic layer deposition. The insulating state shows activated behavior with an energy gap that increases with B. DC source−drain bias measurements reveal a transition between direct tunneling and Fowler-Nordheim tunneling at T = 0.29 and 2 K as bias increases. The results support the existence of a gapped region with a length on the order of ≈100 nm in an ≈440 nm-wide constriction created within a larger area of graphene. This work demonstrates the opportunistic study of properties serendipitously created in a device within a multiplexed array. High-throughput testing via the multiplexing platform increases the opportunity to identify such devices with unusual properties and study them and is compatible with different nanoelectronic devices including 2D materials and semiconductor nanowires, measured from cryogenic to room temperatures. 16

METHODS
These methods concern the primary device on the multiplexing platform. Methods relating to the second device are provided in the Supporting Information. The multiplexer is fabricated in a GaAs high electron mobility transistor (HEMT) in which a 2DEG forms 90 nm below the surface (Cavendish wafer V832). The mean sheet density and electron mobility are n GaAs = 1.53 × 10 11 cm −2 and μ GaAs = 8.17 × 10 5 cm 2 V −1 s −1 , respectively, measured on a separate Hall bar device. A conducting path is defined within the 2DEG from one input to 16 outputs, and individual outputs are chosen by surface addressing gates. Multiplexer layout, fabrication, and operation details are described in ref 16, although the device measured here is different from that in ref 16. Source−drain contacts (Ti/Au) are fabricated on the multiplexer before transferring a monolayer graphene film grown by chemical vapor deposition onto Cu. The graphene is coated with poly(methyl methacrylate) (PMMA), and the Cu is etched using a 1.2% solution of ammonium persulfate etchant (NH 4 ) 2 S 2 O 8 . The PMMA−graphene stack is transferred to the multiplexer after rinsing in a DI water bath and dried in air overnight before baking at 125°C for 30 min and soaking in acetone to remove the PMMA. Graphene channels (10 μm wide) are defined between source−drain contacts (10 μm separation) by photolithography and oxygen reactive ion etching. For the device measured here, Raman measurements show a large 2D to G peak ratio across the entirety of the graphene, consistent with a single layer, shown in Figure S6. The two-terminal differential conductance is measured using 77 Hz AC voltage, V AC = 400 μV (rms), applied to the multiplexer and graphene in series. At B = 10 T, the maximum power dissipated in the graphene and energy loss rate per carrier are estimated to be <0.2 pW and <1.5 × 10 −18 W, respectively, which suggests a carrier temperature of <0.6 K, as discussed in the Supporting Information.
Series resistance, energy gaps estimated using an Arrhenius equation, carrier mobility estimated using the quantum Hall effect, DC source−drain bias across the graphene, effective gating of the graphene due to source−drain bias, DC current−voltage characteristics for negative V DS , energy loss rates and carrier temperature, Raman spectroscopy, and methods for the deliberately defined constriction (PDF)